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ABSTRACT 

We present a technique to extract ultra-deep diffuse-light spectra from the standard 
multi-object spectroscopic observations used to investigate extragalactic globular clus- 
ter (GC) systems. This technique allows a clean extraction of the spectrum of the host 
galaxy diffuse light from the same slitlets as the GC targets. We show the utility of 
the method for investigating the kinematics and stellar populations of galaxies at 
radii much greater than usually probed in longslit studies, at no additional expense 
in terms of telescope time. To demonstrate this technique we present Gemini/GMOS 
spectroscopy of 29 GCs associated with the elliptical galaxy NGC 3923. We com- 
pare the measured stellar population parameters of the GC system with those of the 
spheroid of NGC 3923 at the same projected radii, and find the GCs to have old 
ages > 10 Gyr, [a/Fe]~0.3 and a range of metallicitics running from [Z/H] = -1.8 to 
+0.35. The diffuse light of the galaxy is found to have ages, metallicitics and [a/Fc] 
abundance ratios indistinguishable from those of the red GCs. 

Key words: galaxies: general - galaxies: abundances - galaxies: individual: NGC 
3923 - galaxies: stellar content - globular clusters: general 



1 INTRODUCTION 

In recent years great strides have been made in the spectro- 
scopic study of globular clusters (GCs), and we have begun 
to build up large samples of spectroscopically determined 
ages, metallicities and alpha-element abundances with the 
aim of using these data to infer star formation histories for 
their host galaxies (Kissler-Patig et al. 1998; Kuntschner 
et al. 2002; Puzia et al. 2004, 2005; Strader et al. 2005; Pierce 
et al. 2006a,b; Cenarro et al. 2007; Strader et al. 2007). GCs 
are particularly suitable for this analysis because they can 
be assumed to be simple stellar populations (SSP) where 
all of their stars formed at one epoch from a single cloud 
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of gas with an almost uniform metallicity. This makes them 
much simpler to model than the integrated stellar popula- 
tions of the diffuse light of galaxies, where many star forma- 
tion events at different epochs greatly complicate analysis 
of the spectra. 

The great promise of this approach to simplify the ex- 
amination of the star formation history of galaxies is how- 
ever dependent on one critical assumption: that the stellar 
populations of GCs are representative of the field stars that 
form the bulk of a galaxies' stellar population. This would 
occur naturally if GCs form with some fixed proportion to 
those stars that form in unbound star forming regions dur- 
ing a star formation event and with similar element abun- 
dances. There are several observations that appear to bear 
out this assumption, including observations by Larsen & 
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Richtler (2000) that show that the number of Young Mas- 
sive Clusters (YMCs) correlates with the star formation rate 
per unit area in the galaxies that host them. These YMCs 
are believed to be the low redshift analogues to the star 
clusters that survived for a Hubble time to form the mostly 
old GC populations seen in galaxies today. The almost con- 
stant value observed for the number of GCs normalized to 
the total baryonic mass of a galaxy (McLaughlin 1999) also 
provides strong support for the idea that GCs closely trace 
the major star formation events of a galaxy. 

A more direct way to approach this issue is to examine 
the ensemble parameters of large numbers of GCs and com- 
pare them with those measured for the host galaxies' diffuse 
light. If the GCs do indeed act as a faithful tracer of the dif- 
ferent star formation events in a galaxy's history, then the 
average GC stellar population should closely resemble that 
of the overall galaxy. 

To date little work has been done in comparing the stel- 
lar population parameters measured from line index analy- 
ses for GCs with those determined for the host galaxy dif- 
fuse stellar component. Principally this is because at the 
galactocentric radii where the GCs selected for spectroscopic 
analysis are found, the surface brightness of the underlying 
galaxy halo light is very faint. In Norris et al. (2006) the 
stellar populations along the major and minor axes of the 
isolated SO galaxy NGC 3115 were compared to the stellar 
populations of the GC population of NGC 3115 taken from 
Kuntschner et al. (2002). For this galaxy it was found that 
the stellar population along the minor axis at ~2R e , was 
indistinguishable from that of the more metal rich GCs in 
terms of age (~12 Gyr), [Z/H] (~-0.5 dex) and [a/Fe] (~0.3 
dex). Furthermore, at larger radii along the minor axis the 
stellar population became more metal-poor, implying that if 
it were possible to observe the halo stellar population at even 
larger radii it may begin to resemble the more metal poor 
GCs. These results hint at a strong connection between the 
formation of the spheroid of NGC 3115 and its GC system. 

Ideally we would like to repeat this type of analysis for 
many galaxies to see how universal this behaviour is. How- 
ever the large additional burden of telescope time required 
for the longslit measurements on top of the long exposure 
times required for the MOS (Multi-Object Spectroscopy) 
GC investigations makes this approach impractical. In this 
paper we present a new technique that allows us to combine 
both parts of this type of study into one set of observations 
with no extra cost in terms of telescope time. As an ad- 
ditional benefit, this technique can also be used to extract 
kinematic parameters for the host galaxy at radii well be- 
yond those commonly achieved for standard longslit studies. 

For this study we examine the galaxy NGC 3923 (see 
Table 1 for details), which is a large nearby (D ~ 21Mpc) 
elliptical galaxy with a prominent shell structure (Malin & 
Carter 1983). NGC 3923 is the brightest galaxy of an aver- 
age sized group with a low early-type fraction, perhaps sug- 
gesting that the group is dynamically young (Brough et al. 
2006). 

An HST ACS study of the GC system of NGC 3923 
by Sikkema et al. (2006) confirmed previous observations 
(Zepf et al. 1994, 1995) that the GC system has a bimodal 
colour distribution and an unsually high specific frequency 
(Sat>5) for an early type galaxy in a lower density environ- 
ment. Sikkema et al. (2006) searched for evidence of a GC 



Table 1. NGC 3923 Basic Parameters 



Parameter 



Value 



Right Ascension (J2000) 
Declination (J2000) 
I 
b 

Morphological Type 
Magnitude 
Major Diameter 
Minor Diameter 
Heliocentric Radial Velocity 
J-Band Half-Light Radius 



ll h 51 m 01.8 s 
-28°48'22" 
287.28 
32.22 
E4 

10.8 B mag 
5.9 arcmin 
3.9 arcmin 
1739 ±9 kms- 1 
43.8" f 



Table data from NED: 

http:/ /nedwww. ipac.caltech.edu/ 

Except: 

f 2 Micron All Sky Survey 
www.ipac.caltcch.edu/2mass/ 



sub-population formed during the merger event thought to 
be responsible for the creation of the shell structures. They 
concluded that their photometric data did not support the 
presence of any younger GC population with ages similar 
to that of the shell structures, which they believe to have 
formed 0.8-1.2 Gyr ago. They did however find that the bi- 
modal V-I colour distribution and radial density profiles of 
the blue and red GCs were typical for old GC systems in el- 
lipticals. Sikkema et al. (2007) used the same HST ACS data 
to study the shell structures of NGC 3923 finding them to 
be similar in colour or slightly redder than the galaxy stellar 
light. 

Other studies have examined the stellar population and 
kinematics of the central regions of NGC 3923. Thomas 
et al. (2005) found the nuclear regions to have an age of 
3.3 Gyr, [Z/H] = +0.62 and [a/Fe] = 0.31, while Denicolo 
et al. (2005) using a similar method found values of 2.6 Gyr, 
[Z/H] >0.67 and [a/Fe] = 0.14 respectively. 

This paper is structured as follows. Section 2 provides 
details on the data used to investigate the GC and diffuse 
stellar content of NGC 3923; Section 3 describes the method 
used to extract extra information on the underlying galaxy 
spectrum from a given set of MOS observations; Section 4 
briefly describes kinematic results derived using this tech- 
nique, and provides an example of the line index analysis 
that is possible using this method. In Section 5 we discuss 
the implications of these results, Section 6 provides some 
concluding remarks. 



2 OBSERVATIONS AND DATA REDUCTION 

Whilst we present a specific example using data obtained 
from the Gemini-GMOS instrument, the approach described 
here should be equally applicable to similar MOS intruments 
such as FORS1/2 on the VLT or LRIS/DEIMOS on Keck. 

All observations described here were taken with the 
Gemini South Multi-Object-Spectrograph (CMOS) (Hook 
et al. 2004) as part of Gemini program GS-2004A-Q-9. 

Pre-imaging of NGC 3923 for object selection was un- 
dertaken on 2004 January 19 and consisted of 4x200s expo- 
sures in Sloan g' and 4x 100s in r' and i' for each of 3 fields 
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(one central, one SW and one NE of the galaxy centre). 
For a thorough discussion of the procedure used to reduce 
the pre-imaging and select GC candidates see Forbes et al. 
(2004) and Bridges et al. (2006) . A full examination of the 
photometric properties of the GC system of NGC 3923 and 
other galaxies investigated in this project will be presented 
in a forthcoming paper Faifer et al. (in prep). 

GMOS Nod and Shuffle masks were produced for each 
of the three fields in addition to one MOS mask for the cen- 
tral pointing (see Fig. 1). Since it is the MOS mask that 
provides the data described here, we shall leave any further 
discussion of the Nod and Shuffle masks and the kinematic 
investigation of NGC 3923 derived from them to a later pa- 
per (Norris et al. in prep.). 

The MOS mask consisted of 37 slitlets of width 1 arc- 
sec by a minimum length of 4 arcsecs. The mask was ex- 
posed using the B600_G5303 grism for 8 x 1800s at a cen- 
tral wavelength of 500nm and 8 x 1800s at a central wave- 
length of 505mn (to cover the CCD chip gaps), yielding 
8 hours of on-source integration. The seeing ranged from 
0.6 - 0.9 arcsec during the observations. Bias frames, flat 
fields and copper-argon (CuAr) arc spectra were observed 
throughout the observations as part of the Gemini baseline 
calibrations. The MOS spectra produced typically cover the 
wavelength range 3900-5500 A although because the wave- 
length range depends on slit position, some spectra start at 
~3500A while others end at ~7200A. 

Data reduction of the MOS spectra through to the point 
of producing rectified and wavelength calibrated 2D spec- 
tra was accomplished utilising the Gemini/CMOS packages 
in IRAF as described in Bridges et al. (2006). From the 
CuAr arcs, wavelength calibrations with residuals ~0.lA 
were achieved. The extraction of the target GC spectra 
from each 2D spectrum was undertaken using the APALL 
task in the APEXTRACT package using an optimal (vari- 
ance weighted extraction) . Flux calibration was achieved us- 
ing GMOS longslit observations of the flux standard star 
LTT3864 made during the same semester with identical ob- 
servational set-up. 

The extraction of the NGC 3923 diffuse light spectra 
from the MOS slitlets was undertaken using a custom 1DL 
script implementing the algorithm described below. 



3 METHOD 

Modern multi-object spectroscopic instruments such as 
GMOS are capable of taking spectra of around 40 objects 
simultaneously over fields of view of around 6x6 arcmin 2 . 
Each MOS slitlet is typically around 1 arcsec wide by a few 
(4-10) arcsec long. Generally the GCs are unresolved from 
the ground producing spectra smeared over ~1 arcsec (in 
the spatial direction) for typical seeing values. The rest of 
the slitlet collects background (sky + galaxy) photons used 
to background subtract the GC spectrum. 

For the typical field of view of the current generation 
of MOS instruments and the distance and size of the tar- 
get galaxies being investigated by these studies, it can be 
seen that significant amounts of target galaxy flux are con- 
tained in the MOS slitlets being used to study the GCs (See 
Fig. 1) . The flux incident on each slitlet can then be thought 
of as a sum of contributions from the target GC, the actual 
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Figure 1. NGC 3923 GMOS i' band image. Circles show the posi- 
tion of photometrically chosen targets later confirmed by spectro- 
copy to be GCs, stars show targets found to be Milky Way stars 
and squares are targets found to be background galaxies/QSOs. 
The problem of placing a slitlet at large enough galactocentric 
radii to produce an uncontaminated sky spectrum when using 
modern MOS instruments with fields of view around 6x6 arcmin 2 
is apparent. Also visible are the shell structures of NGC 3923 on 
the major axis. 



sky background (both atmospheric and extra-galactic) and a 
contribution from the diffuse light of the target galaxy. Fig. 2 
shows an example of a typical 2D MOS spectrum, showing 
the regions containing the different contributions to the flux. 
If the instrument has a sufficiently large field of view, it is 
possible to place a slitlet at large galactocentric radii where 
the spectrum measured would be essentially a pure sky spec- 
trum, uncontaminated by galaxy flux. This sky spectrum 
could then be subtracted from the background regions of 
each of the GC slitlets to produce a sky-subtracted galaxy 
diffuse light spectrum for each of the individual slitlets. We 
have successfully utilised this approach in a study similar 
to this one involving the GC system of NGC 524, finding 
that slitlets located beyond 5R e provide an adequate sky 
spectrum. 

In practice the FOV of the current generation of MOS 
instruments is often too small to allow a slitlet to be located 
at large enough galactocentric radii to produce a "pure" 
sky spectrum of the accuracy required (at least for stellar 
population studies). 

Here we present an algorithm that allows the calcula- 
tion of the sky spectrum even in cases with small FOV, by 
combining the spectroscopic dataset with wider field imag- 
ing data. 



3.1 Algorithm 

• The spectra are reduced as normal to the point where 
the individual spectra are cut out, rectified and wavelength 
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Figure 2. Example of a typical 2d MOS spectrum, showing tar- 
get spectrum (in this case a GC spectrum) and sky regions typi- 
cally used to sky subtract the target spectrum. For typical MOS 
GC studies these sky regions can contain useful amounts of host 
galaxy diffuse light. 

calibrated. See Bridges et al. (2006) for more details of the 
reduction used for the data presented here. 

• A sky + galaxy spectrum is produced for each slitlet. 
The spectra are produced using the usual method for cal- 
culating a sky spectrum from a 2D slitlet (i.e the spectrum 
used to sky subtract a 2D frame). Fig. 3 displays three such 
sky + galaxy spectra for a range of galactocentric radii; as 
is expected those spectra taken from slits which lie closer to 
the galaxy centre contain more flux. 

• Each spectrum is rebinned onto a uniform wavelength 
scale. 

• The wider field pre-imaging is combined to produce a 
master wide field image of the target galaxy. In this case we 
make use of our g' band imaging as this best matches the 
wavelength coverage of our spectra, though we note that if 
our i' band imaging is used instead we obtain line indices 
which differ by less than 0.1 A. We caution however that 
care should be taken to ensure that observed optical light 
profile being used should match that of the spectral region 
being examined, i.e. that the use of an IR image to set the 
scaling for optical spectra is probably unwise because of the 
different photometric profiles displayed by galaxies in the 
different wavelength regimes. 

• Using the master image of the field, the local back- 
ground in an annulus around each of the target GCs is cal- 
culated. 

• The actual sky background number of counts in the 
image is measured for a region located at large galactocen- 
tric radii. Experience with the data for NGC 3923 showed 
that this value was always very similar to the local back- 
ground measured for the GC with the largest galactocentric 
distance. 

• For each wavelength pixel of the rebinned spectra the 
number of counts is plotted against the local background 
around the corresponding object determined from the wide 
field image. (See Fig. 3 for the pixels corresponding to 
5100 A). As can be seen in the lower panel of Fig. 3 the 
correlation of local background with number of counts in 
the spectra is reasonably good. This correlation can then 
be fit with a simple straight line and extrapolated to the 
value calculated for the sky background at large radius from 
the wide field image. In the example in Fig. 3 the value of 
the sky brightness measured at large radii to represent the 
uncontaminated sky value was around 1000 counts in the 
g' band, leading to an estimate of the "pure" sky spectrum 
having around 32 counts in the pixel located at 5100 A . 

• When repeated for each wavelength pixel, this proce- 
dure results in a spectrum that is a good approximation of a 



"pure" sky spectrum. Note that the produced sky spectrum 
is only accurate in the regions where a significant number 
of the spectra overlap in wavelength coverage. In the case 
of our GMOS data this restricts our analysis of the galaxy 
spectra to the wavelength range 4400 - 5500 A where the 
majority of the 37 input spectra have coverage. 

• Each of the 37 rebinned sky + galaxy spectra are then 
sky subtracted using the "pure" sky spectrum, leaving 
37 spectra that are essentially pure galaxy spectra. The 
spectra measured from individual exposures can then be 
combined in the usual manner; examples of 3 of the spectra 
produced by combining 16 individual exposures can be seen 
in Fig. 4. 

In cases where wider field imaging is not available it is pos- 
sible to obtain similar results using an iterative procedure. 
In this case an initial guess for the sky brightness assumed 
to represent pure night sky is made (by extrapolating the 
observed galaxy profile for instance) and the entire proce- 
dure from steps 6 to 8 can then be iterated with the assumed 
night sky brightness being varied until sky line residuals in 
the final diffuse light spectra are minimised. 



3.2 Kinematics 

GC radial velocities were measured using the FXCOR task 
in the RV package of IRAF. We did not observe radial ve- 
locity standard stars as part of this investigation due to the 
extra overhead in observing time this requires. For templates 
we therefore use spectra from the simple stellar population 
library of Vazdekis (1999). These templates span a range of 
age of 1-18 Gyr and a metallicity range from [Fe/H] = +0.2 
to -0.7. To this sample we add 6 stellar spectra from the 
Jones (1997) library to extend the metallicitiy range cov- 
ered by the templates to [Fe/H] < -1.5. This wide range 
of age and metallicity helps to minimise the effects of tem- 
plate mismatch on the measured velocities. We use a 3 sigma 
clipped mean of the velocities derived from these templates 
as the final velocity for each GC. The errors are estimated 
from the mean of the errors measured by FXCOR for those 
velocities not removed by the clipping procedure. 

Objects with velocities in the range 1200-2400kms _1 are 
assumed to be associated with NGC 3923 and were classified 
as GC. We find that 29 out of 37 (78%) objects observed 
with our MOS mask have velocities within this range. Of the 
remaining objects, 2 were z~0.3 galaxies, one was a z~1.3 
quasar and 5 were Milky Way stars. 

The recessional velocities of the NGC 3923 diffuse light 
spectra were measured using the publicly available IDL im- 
plementation of the pPXF method (Cappellari & Emsellem 
2004). This method was chosen as it allows a simultaneous 
measurement of the radial velocity and velocity dispersion 
from the spectrum, both of which are important in studies of 
the diffuse light of galaxies. In this instance the spectra were 
fit in pixel space over the wavelength range 4800-5400 A us- 
ing the same templates as the GCs. The final velocity and 
velocity dispersion for each slitlet is taken to be the mean 
of the velocities/velocity dispersions determined from 100 
Monte-Carlo resimulations of the input spectra with added 
photon noise. The error in the measured velocity is taken to 
be the dispersion in the 100 best fit velocities. As a check 
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Figure 3. Upper Panel: Upper three spectra are the "sky + 
galaxy" spectra extracted from slits covering the range of galac- 
tocentric radii covered in this study. The fourth spectrum (black 
line) shows the extrapolated "pure" sky spectrum using the algo- 
rithm described in the text. Mid Panel: This spectrum is the result 
of subtracting the extrapolated sky spectrum from the uppermost 
"sky + galaxy" spectrum. Redshiftcd (z=0.006) absorption lines 
due to H/3, Mgb, Fe5270 and Fe5335 are clearly visible. Lower 
Panel displays a typical result of the extrapolation process used 
to determine the actual pure sky spectrum in this case for the 
pixel located at 5100 A. The position of the 3 points from panel 
A is also displayed here for instructive purposes, as are the mea- 
sured number of counts and statistical error at this pixel (5100 A ) 
for the other 34 slitlets (black circles) in this mask. 
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Figure 4. Representative diffuse light spectra for NGC 3923 mea- 
sured using our new technique. These spectra are the results of 
coadding 16 individual exposures, they have been smoothed to 
the Lick/IDS resolution. Redshifted absorption lines due to sev- 
eral species are present. 



of consistency we remeasured the recessional velocity of the 
diffuse light spectra using the FXCOR task in 1RAF and 
found no significant differences between the two methods. 

3.3 Line Indices 

By binning all of our 37 diffuse light spectra together we have 
produced one spectrum of sufficient signal-to-noise (S/N > 
60 per A) to measure accurate absorption line indices. Since 
NGC 3923 shows evidence of slight rotation (Fig. 5) it is 
necessary to remove the effects of rotation on the spectra 
by correcting the spectra to some common velocity before 
coadding them. In this instance the velocity of each of the 
37 spectra was measured as described previously, then using 
the DOPCOR task in 1RAF each of the spectra was de- 
redshifted to zero redshift before combining. 

For ease of comparison with other studies we have cho- 
sen to use the Lick/IDS system when examining the stellar 
population present in this master spectrum and in our GC 
sample. There are several steps that must be taken to ensure 
that our measured indices are securely on the Lick/IDS sys- 
tem (Worthey & Ottaviani 1997; Trager et al. 1998). The 
first is to convolve our spectra with a wavelength depen- 
dent Gaussian kernel to reproduce the variable Lick/IDS 
resolution of 9-11 A . We then measure the Lick/IDS indices 
for our fluxed spectra using the wavelength definitions from 
Worthey & Ottaviani (1997), Trager et al. (1998) and the 
INDEXF index measuring code of N Cardiel 1 . Uncertainties 
on the measured indices are a combination of errors intro- 
duced by photon (both target and sky) noise and errors in 
the measured redshift of the spectra (essentially negligible 
in all cases). 

The line indices for the spheroid spectra must also be 
corrected for the smearing effects of the line-of-sight veloc- 
ity distribution (LOSVD). This effect is entirely negligible 
for GC spectra with velocity dispersions of the order a few 
kms -1 , but cannot be ignored for galaxies with velocity dis- 
persions over ~100 kms -1 . As our measured diffuse light 
spectrum has a velocity dispersion of around 200 kms - this 
step is important. We correct for this effect using the proce- 
dure of Kuntschner (2004) which can also be used to correct 
line indices for the additional influences of the higher order 
LOSVD terms I13 and I14. Kuntschner (2004) provides a new 
parameterisation for the LOSVD line strength corrections 
derived by determining the difference in measured Lick/IDS 
line strengths measured from template spectra before and 
after convolving with various choices of a, I13 and I14. For 
NGC 3923 the correction for li3 and I14 is not implemented 
due to the effect of the correction given the measured values 
of a, I13 and I14 being minimal relative to the other errors. 

The final step required to firmly fix our measurements 
to the Lick/IDS system is to measure standard stars from 
the Lick library using our observational set-up and com- 
pute offsets between the two. Although we did not observe 
Lick standard stars directly during this project, due to the 
prohibitive cost in observing time, data collected as part of 
GS-2003B-Q-63 to calibrate GMOS-S onto the Lick system 
were kindly made available by Bryan Miller. We have used a 

1 Available at: 

http: / / www. ucm.es /info/ Astrof / software / indexf / indexf . html 
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set of 28 stars from this sample, taken with an almost identi- 
cal set-up to our MOS observations and used them to derive 
offsets for GMOS-S to the Lick system. We have applied 
these offsets to our measured line indices presented in Table 
2. We note that while useful for comparison between obser- 
vations taken using different instruments, or for comparing 
measured indices to models based on the Lick/IDS system 
e.g. Thomas et al. (2003, 2004), this step is not essential 
for a simple first order differential comparison between two 
sets of observations made with the same instrument, such 
as that described here. 

Using the procedure outlined (neglecting the higher or- 
der LOSVD corrections) we have measured Lick/IDS line 
strength indices for all 29 of our confirmed GCs, as well as 
our single coadded spheroid diffuse light spectrum. This sin- 
gle diffuse light spectrum is located at a luminosity weighted 
distance of 3 R e (throughout we use R e =43.8 arcsec mea- 
sured in the J-band (Jarrett et al. 2003), though we note 
that V, R and I band estimates from Bender et al. (1988) 
are essentially identical). The S/N of our GC spectra range 
from ~ 11-56 per A measured in the region 5100-5150 A, 
giving errors in the H/3 index of 0.12-0.72 A . 



4 RESULTS 

Table 2 presents velocities and Lick/IDS line strength in- 
dices for all 29 of our confirmed GCs, as well as the single 
coadded diffuse light spectrum, and a series of GC composite 
indices binned by g'-i' colour. In addition to examining the 
GC spectra we have extracted 37 diffuse light spectra from 
our MOS exposure (see Fig. 4 and Table 4) with signal-to- 
noise ranging from 5 to 12, which is sufficient to measure 
velocities in all cases and velocity dispersions in most cases. 



4.1 Kinematics 

A more detailed discussion of the kinematics of NGC 3923 
and its GC system will be presented in Norris et al. (in prep), 
where we include GC velocities from additional masks at 
larger radii obtained using the Nod and Shuffle technique 
(Cuillandre et al. 1994). 

In this section we limit ourselves to examining the 
spheroid of NGC 3923 for rotation and in Section 4.2 to com- 
paring the velocities measured for the spheroid of NGC 3923 
and those measured for the GCs examined in the same slits. 

Here we present the results of our search for evidence of 
rotation in the diffuse stellar halo of NGC 3923 (see Fig. 5) 
by carrying out a non-linear least squares fit to the equation: 

V(0) = V rot sin(6» - 6» ) + V 

where V(0) is the velocity measured for the diffuse 
light in each slitlet, V ro t is the amplitude of the projected 
rotation velocity, Vo the systematic velocity of NGC 3923, 
9 the azimuthal angle of each slitlet relative to the galaxy 
major axis and 6q is the position angle of the line of 
nodes. This approach determines the best-fitting simple fiat 
rotation curve (see Zepf et al. (2000) for details). 

We find evidence for rotation along the major axis of 
NGC 3923, with the best fitting amplitude being 31 ± 13 



kms^ 1 along a position angle of 290 ± 19° where the ma- 
jor axis is approximately at PA = 315°. This result is in 
good agreement with that of Koprolin & Zollinger (2000) 
who found a small (53 ± 13 kms^ 1 ) but non negligible ma- 
jor axis rotation for NGC 3923. It is however different to 
the results of Carter et al. (1998) who found that the in- 
ner 25 arcsec of NGC 3923 showed no rotation on its major 
axis but had minor axis rotation of amplitude ~ 20 kms -1 . 
The difference in conclusions between the different studies 
may be attributed to the different radial ranges over which 
they are sensitive; our method provides a measurement of 
the rotation of the galaxy at a radius at least a factor of 
two greater than the largest radial point available to either 
Carter et al. (1998) or Koprolin & Zeilinger (2000). 

These measurements of the kinematics of NGC 3923 
cover a range in effective radii of between 2 and 4 R e , demon- 
strating the great power which this approach has for illumi- 
nating the kinematics of the outer regions of galaxies. 

4.2 Reliability of Velocities 

One concern with our procedure for extracting diffuse host 
galaxy spheroid light from the same slitlets as the target GC 
spectra is that the spheroid spectra could be contaminated 
by flux from the target GCs. In general this should not be of 
major concern due to the careful selection of the sky regions 
of the 2D spectrum. As a test of this assumption in Fig. 6, 
we plot the measured diffuse light velocity against the target 
velocity (GC, Milky Way star, background galaxy or QSO) 
from the same slitlet. In Fig. 6 the filled circles are target 
objects classified by their velocities as GCs, filled stars are 
objects classified as Milky Way stars, the three background 
objects (two z~0.3 and a z~1.3 QSO) are omitted from the 
plot for clarity, although they are included in the fit which 
produces the dashed best fit trend line (which has slope 
equal to zero within the errors) . As can be seen clearly, there 
is no correlation between the velocities measured for the 
target spectra and those measured for the diffuse light; in 
particular even for the brightest target objects (generally the 
MW stars) the measured velocity for the spheroid of NGC 
3923 is still consistent with the recessional velocity of NGC 
3923. This indicates that contamination of the diffuse light 
spectra by the target object spectra is negligible. 

4.3 Stellar Populations 

In this section we examine the stellar populations of the GCs 
and the diffuse light of NGC 3923. To begin we shall discuss 
some of the qualitative properties of the stellar populations 
through examination of age-metallicity diagnostic plots. We 
shall then examine the data more quantitatively by present- 
ing the ages, metallicities and a-element abundances derived 
through x 2 fitting of the data to stellar population models. 

4-3.1 Index-Index Plots 

In the left panel of Fig. 7 we present an Mgb vs. <Fe> 
plot 2 for the 29 GCs as well as the NGC 3923 diffuse light 
spectrum measured in this study. Model predictions from 

2 <Fc> = (Fc5270 + Fc5335)/2. (Gonzalez 1993) 
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Figure 5. Velocity vs. Azimuthal angle for NGC 3923 spectra 
extracted from MOS slitlets (circles). The solid line is the best fit 
flat rotation curve indicating rotation amplitude of 31 kms~ 1 ± 
13 with position angle 290° ± 19 (north through east). 
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Figure 6. Measured velocities from NGC 3923 diffuse light spec- 
tra vs. the measured velocities of the targets observed in the same 
slits. Stars indicate targets classified as Milky Way stars, filled cir- 
cles those targets classified as GCs belonging to the NGC 3923 
system. Two z~0.3 galaxies and a z~1.3 quasar arc excluded for 
clarity. The dotted line shows a linear fit to the data, indicating 
no correlation between the velocities measured for the targets and 
the galaxy diffuse light. 

Thomas et al. (2003, 2004) are overplotted for [a/Fe] = 0.0, 
0.3 and 0.5, ages 1-12 Gyr and metallicity running from - 
2.25 (bottom left) to +0.67 (top right). This choice of indices 
produces a figure which is insensitive to age but can be used 
to constrain the metallicity and a-element abundance ratios 
(Kuntschner et al. 2002). As can be seen the GCs (light grey 
circles) have metallicities ranging from -1.7 to +0.5 and ap- 
pear to have a wide range of a-element abundance ratios. 



When the GC data is binned by colour (coloured circles) 
the behaviour of the binned data becomes more systematic, 
with the points showing a clear preference for a constant a 
element abundance ratio of ~0.3, even as the metallicity of 
the binned data decreases from around 0.0 to -1.5. These 
bins cover the range of colour displayed by the majority of 
extragalactic GCs, and have mean g' - i' colour of 0.82, 0.99, 
1.09 and 1.19 (See Table 2). The position of the mean colour 
of each bin relative to the peaks of the blue and red GC 
distributions can be seen in Fig. 8. In particular the mean 
colour of the third bin can be seen to correspond closely to 
that of the peak of the red GC distribution. If all the GC 
data are combined to produce a composite GC data point 
(white triangle); the combined GC spectrum has a metal- 
licity of ~ -0.8 and an abundance ratio consistent with 0.3 
within the (significant) errors. It is interesting that both 
the central regions of NGC 3923 measured by Trager et al. 
(1998) and Beuing et al. (2002) (grey triangle and square) 
and the NGC 3923 diffuse light spectrum measured in this 
study (white star) also lie along the same [a/Fe] = 0.3 track. 
However the fact that the data point from Denicolo et al. 
(2005) (black triangle) for the central R e /8 of NGC 3923 is 
offset relative to the other two observations illustrates the 
magnitude of possible systematic errors in this type of study. 
The difference in metallicity between the inner nuclear data 
and our large radii diffuse stellar data point can be under- 
stood in terms of a radial metallicity gradient in the galaxy. 
A further point of note is the close agreement between the 
mean of the GC spectrum and the NGC 3923 diffuse light 
spectrum which was measured at the same projected radii. 

The right panel of Fig. 7 presents an [MgFe]' vs H/3 age- 
metallicity plot 3 with symbols as previously defined. Here it 
is necessary to choose a fixed value of [a/Fe] for the grids 
and we have chosen 0.3 as this value is close to that found 
for each of the populations described in the previous sec- 
tion. It can be seen that three GCs appear to lie along grid 
lines consistent with young or intermediate ages, these ob- 
jects are however generally among the lowest S/N objects in 
our sample. When other age-metallicity diagnostic plots (H7 
and US vs. [MgFe]') are considered as well as the fact that 
similar numbers of objects lie below the grid it appears that 
these objects may just be statistical outliers in a distribu- 
tion centred on uniformly older ages. We shall return to this 
point when examining the results of the \ 2 fitting approach. 
The binned GC data again behaves in a systematic manner, 
in which all the data points lie along a constant age line of 
around 12 Gyr, with metallicities consistent with those seen 
in the previous plot. The mean of the GC system as a whole 
and the NGC 3923 diffuse light point can also be seen to be 
very consistent in terms of age and metallicity. 

The close agreement in apparent age, metallicity 
and [a/Fe] between the NGC 3923 diffuse light spectrum 
measured at 3 R e and the GC system (especially the 
red GCs) at the same projected radii is tantalising. Such 
close agreement hints at a deep connection between the 
two populations exactly as would be expected if the two 
populations formed coevally from the same material. This 
interpretation is however dependent on the assumption that 



3 [MgFe]' = y/Mgb x (0.72 x Fe5270 + 0.28 x Fe5335 (Thomas 
et al. 2003) 
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Figure 7. Left Panel. Comparison of the [o/Fe] ratios of the diffuse light and GC population of NGC 3923 through the use of an Mg6 
vs. <Fe> diagram. Light grey circles are the 29 GCs measured here (error bars omitted for clarity). The coloured circles are the result 
of binning the 29 GC data points into 4 bins of 7/8 GCs per bin by g'-i' colour, the error bars show the la scatter in each bin. The 
white triangle is the result of binning all 29 GC spectra into one bin, again the error bar represents the la scatter in the bin. The grey 
triangle, square and black triangle are the central NGC 3923 values from Trager et al. (1998), Beuing et al. (2002) and Denicolo et al. 
(2005) respectively. The white star shows the values measured from the single high S/N coadded spectrum of the diffuse light of NGC 
3923. Overplotted are models by Thomas et al. (2003, 2004) with abundance ratios [a/Fe] = 0.0, 0.3, 0.5, the models have ages 3-12 Gyr 
and metallicity [Z/H] = -2.25, -1.35, -0.33, 0.0, +0.35 and +0.67. Right Panel. Age-Metallicity diagnostic plot. Symbols as in left panel. 
Models are overplotted for a fixed value of [ct/Fe] = 0.3 



the GCs examined here are representative of the full GC 
population. We believe that our sample is a fair match to 
that present in NGC 3923 over the range of projected radii 
studied, though only observation of a significantly larger 
sample of GCs carefully chosen to match the full sample in 
terms of colour and luminosity would allow us to determine 
this for certain. 



4-3.2 Derived Stellar Population Parameters 

We derive ages, metallicities and Q-element abundances for 
the stellar populations being studied here using the multi- 
index fitting method of Proctor & Sansom (2002) and Proc- 
tor et al. (2004). This technique compares the measured Lick 
indices with SSP models, in this case the SSP models of 
Thomas et al. (2003, 2004) (TMK04) which to date remain 
the only SSP models in which the effects of a-element abun- 
dance ratios on Balmer line indices at low metallicity are 
included. For a more thorough description of the different 



treatments of Q-element abundance ratios in modern SSP 
models see Mendel et al. (2007). 

The procedure used to determine the best fitting SSP 
model is as follows. The total set of measured Lick/IDS in- 
dices is compared to the TMK04 SSPs and a minimum \ 2 
fit obtained; simultaneously a set of \ 2 minimization fits are 
determined with each of the indices omitted. The lowest to- 
tal \ 2 fit from this set is then chosen, the necessary index 
removed and the process repeated until a stable fit is ar- 
rived at where no highly aberrant (3cr) indices remain. In 
this instance all GCs displayed some indices which were suf- 
ficiently aberrant to be excluded from the fit, though all fits 
utilised a minimum of 9 indices, most GCs are fit with at 
least 13 indices. Errors in derived parameters are determined 
via monte-carlo re-simulations of the input index measure- 
ments with their measured errors. The results of this proce- 
dure are displayed in Table 3. 

In common with previous studies we find that the 
molecular bands Mgi and Mg2 appear depressed, an ob- 
servation generally explained as being due to problems with 
flux calibration. We have tested this explanation by exam- 



© 2007 RAS, MNRAS 000, 1-15 



Globular Clusters and the Spheroid of NGC 3923 9 



ining the effect of using the flux calibration curves found 
from repeat observations of the same standard star (in this 
case LTT4364) observed for a different program using the 
same observational set up used here. We found that even 
when using curves derived from the same standard star the 
measured Mgi and Mg2 indices of our GCs after flux cali- 
bration could vary systematically by as much as 0.01 mag, 
which is of the same order as the deviations observed in our 
fitting procedure. Other indices were much less affected by 
the change in flux calibration, presumably due to the much 
narrower wavelength coverage of these indices. Accordingly 
Mgi and Mg2, as well as Ca4227 which was similarly de- 
pressed, were therefore removed from our fitting procedure. 
The CN indices tended to be enhanced relative to the mod- 
els, in agreement with previous findings (Beasley et al. 2004; 
Pierce et al. 2006a) for GCs. The Ca4668 index also tended 
to be unreliable and was therefore excluded from the fits. 

Our quantitative x 2 fitting procedure produces results 
in good agreement with the more qualitative discussion in 
the previous section. We find the GCs examined to have ages 
consistent at the one sigma level with old ages ( > 10 Gyr) 
in all but one case. The GCs exhibit a spread of metallici- 
ties from [Z/H] = -1.8 to +0.35 and a-element enhancement 
ratios generally consistent with ~0.3 dex, although several 
objects display [a/Fe] of below zero or greater than 0.5, at 
the limits of the models. From our sample of 29 GCs only 
one object (ID: 197) displays an age which is inconsistent 
with a mean age greater than 10 Gyr, however this result is 
not significant being only a 1.5a deviation. A more thorough 
analysis of the distribution of stellar population parameters 
is somewhat limited by the small sample size, however to 
first order we can state that the distribution of ages is consis- 
tent with a single value of around 12Gyr plus measurement 
errors. The distribution of metallicity however is inconsis- 
tent at the greater than 99.9% significance with having a 
single value plus measurement errors. We find that the dis- 
tribution of metallicity measured is consistent with a mean 
of -0.65 dex with an intrinsic scatter of 0.5 dex plus measure- 
ment errors, assuming a gaussian distribution of metallici- 
ties. Similarly the distribution of [a/Fe] is also inconsistent 
with a single value plus measurement errors at the greater 
than 99.8% significance level, in this case we estimate an in- 
trinsic scatter of 0.12 dex plus measurement errrors around 
a mean of 0.31 dex. 

We have repeated this analysis for our binned GC data 
and for the line indices measured for the diffuse light halo 
of NGC 3923. We find that the binned GC data again be- 
haves in a manner consistent with the more qualitative ex- 
amination of the index-index plots. All bins have ages con- 
sistent with being older than 10 Gyr, all four bins display 
a smoothly increasing metallicity with colour, ranging from 
[Z/H] = -1.375 to -0.125, and all display [a/Fe] consistent 
with 0.3. The considerably larger errors present in the mea- 
sured parameters for the reddest colour bin, are a result 
of this bin being composed almost entirely of lower lumi- 
nosity and hence lower S/N objects. Apparently through 
chance selection this bin contains GCs of considerably lower 
lumnosity than the average for GCs examined spectroscop- 
ically here (see Fig. 8). This results in a mean g' magnitude 
fully half a magnitude below any of the other bins; together 
with the effect of error weighting the indices when combin- 
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Figure 8. Top Panel. Light grey histogram shows the distribu- 
tion of g'-i' colours for photometrically selected candidate GCs 
around NGC 3923 from our GMOS pre-imaging. Dark histogram 
(multiplied by a factor of 10 for clarity) shows the distribution of 
g'-i' colours for our 29 spectroscopically confirmed GCs. Bottom 
Panel, g'-i' vs g' colour-magnitude diagram for our photomet- 
ric (light grey circles) and spectroscopically (dark grey circles) 
confirmed GCs. The coloured circles indicate the error weighted 
mean values for each bin, the dashed lines indicate the arbitrary 
choice of bin ranges. 

ing them this explains why this bin is significantly lower S/N 
than the others. 

Our "Super GC" bin, containing the error weighted in- 
dices of all 29 GCs unsurprisingly predicts ages, metallicities 
and a-element abundances in good agreement with those of 
the colour binned GCs. We derive an old age for the GC 
population as a whole, a metallicity [Z/H] = -0.8 and an 
a-element abundance ratio which is consistent with 0.3. 

Turning our attention to the measured parameters for 
the diffuse light of NGC 3923 at ~3R e we again find old 
ages, metallicity [Z/H] of -0.325 and [a/Fe] of 0.34. These 
results are entirely consistent with that found for the red 
bin of GCs (see Fig. 9), those GCs displaying g'-i' colours 
in the range 1.02 to 1.14. 

4-3.3 Stellar Population Parameter Corellations 

In Fig. 9 we search for correlations between the ages, metal- 
licities and alpha-element enhancements measured for the 
various stellar populations using the \ 2 fitting procedure 
outlined previously. 

One interesting observation is that we see no evidence 
for the observed trend of decreasing [a/Fe] with increasing 
metallicity seen in the GC systems of some galaxies (e.g. 
NGC 3379 (Pierce et al. 2006a), NGC 4649 (Pierce et al. 
2006b)), but not in others (NGC 1407 (Cenarro et al. 2007), 
VCC1087 (Beasley et al. 2006)). 

Another interesting feature visible in our higher S/N 
data (filled circles) is that the Thomas et al. (2003, 2004) 
models tend to introduce a slight trend of increasing age 
with increasing metallicity. This trend, which has been pre- 
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viously observed in the NGC 1407 system (Cenarro et al. 
2007) and the MW GC system (Cenarro et al. 2007; Mendel 
et al. 2007) appears to be a peculiarity of the Thomas et al. 
(2003, 2004) models; Mendel et al. (2007) found that the Lee 
& Worthey (Lee & Worthey 2005) and Vazdekis (Vazdekis 
2007) SSP models do not display this trend. However as the 
age determinations at the lowest metallicities are already 
highly uncertain due to the possible effects of horizontal 
branch morphology, this effect is not significant and does 
not affect our conclusions. 



5 DISCUSSION 

In summary, we have obtained spectra of 29 GCs associ- 
ated with NGC 3923, as well as 37 spectra of the diffuse 
light at radii of 2-4 R e . By combining these 37 diffuse light 
spectra we have produced a single spectrum with a luminos- 
ity weighted distance of 3 R e of sufficient S/N to determine 
an age, metallicity and a-element abundance ratio by com- 
parison with Lick line strength indices for Simple Stellar 
Population models. 

We find evidence for slight rotation of amplitude 31±13 
kms" 1 in the spheroid of NGC 3923 between 2 and 4 R e 
approximately along the major axis of the galaxy. 

All 29 GCs examined in this study are consistent with 
being old ( > 10 Gyr) and none display any persuasive evi- 
dence for being young enough to be associated with the shell 
forming event. Our results therefore support the conclusions 
of Sikkema et al. (2006) who, using HST ACS photometry 
for the GC system of NGC 3923, found no evidence of a 
young GC subpopulation. 

The GCs examined display a wide range of metallicity 
from [Z/H] = -1.8 to +0.35. In contrast to other studies 
(Puzia et al. 2005; Pierce et al. 2006a, b) wc find no evidence 
for a significant subpopulation of GCs with extremely high 
([a/Fe]>0.5) a-element abundance ratios. Our data also do 
not show any clear tendency for a-element abundance to in- 
crease with decreasing metallicity as had been found in these 
studies. It must be noted however that a-element abundance 
ratios are not well constrained at low metallicity and that 
the results of these previous studies may be consistent within 
the errors with a universal value of [a/Fe] = 0.3. In agree- 
ment with both these studies and with measurements of the 
Milky Way GC population we find that the majority of our 
GCs with accurate abundance ratios have [a/Fe] consistent 
with 0.3. 

The stellar population of the diffuse light of NGC 3923 
between ~ 2 and 4R e is entirely consistent in age, metallicity 
and a-element abundance ratio with the red GC population. 
It also agrees within the respective errors with the mean 
of the GC population examined spectroscopically; however, 
there is a slight difference in metallicity for the mean of 
the total GC population and the stellar population of the 
diffuse light of the galaxy. With the current sample size, it 
is impossible to determine if this offset is real or merely the 
result of a statistical fluctuation. 

This similarity between the stellar populations of the 
spheroid of NGC 3923 at large radii and the redder GCs is in 
good agreement with previous findings, both spectroscopic 
(in NGC 3115; Norris et al. (2006)) and photometric (in 



NGC 4649; Forbes et al. (2004)), and supports the view 
that the two populations formed coevally. 

The fact that the stellar populations of the GCs studied 
here do not match those of the inner regions of NGC 3923 
is most likely due to the fact that GC systems are observed 
to display variation in the mix of red and blue populations 
with radii. Therefore the GC populations which make up 
the region between 2-4R e are not likely to be the same as 
those that inhabit the inner regions of the galaxy. Addition- 
ally GCs are thought to be vulnerable to disruption by tidal 
effects within the inner regions of galaxies meaning that pro- 
jection effects become a more serious issue at smaller galac- 
tic radii, with GCs at large radii projected onto the inner 
regions of the galaxy becoming more important. 

By comparing our stellar population parameters for the 
spheroid of NGC 3923 at large galactocentric radii with 
those of the inner regions of NGC 3923 measured by Thomas 
et al. (2005) (using the index measurements of the inner 
Re/10 by Beuing et al. (2002)) and Denicolo et al. (2005) 
(using index measurements of the inner R e /8), we find ev- 
idence for metallicity and age gradients in the spheroid of 
NGC 3923 but no evidence of an [a/Fe] abundance ratio 
gradient. Specifically Thomas et al. (2005) find a luminos- 
ity weighted age for the central region of the galaxy of 3.3 
Gyr, a metallicity of [Z/H] = 0.62 and [a/Fe] = 0.31. This 
is in good agreement with the conclusions of Denicolo et al. 
(2005) who find an age of 2.6 Gyr, metallicity of [Z/H] > 
0.67 and a slightly lower value of [a/Fe] of around 0.14. This 
close agreement is to be expected as both studies make use 
of the same TMK04 models to derive their stellar population 
parameters. We see therefore that the outer regions of NGC 
3923 are considerably older, less enriched in metals, but have 
a similar a-element abundance ratio as the inner regions of 
the galaxy. The first two points are not unexpected consid- 
ering that the values for Thomas et al. (2005) and Denicolo 
et al. (2005) are provided for apertures of < R e /8, and our 
spectroscopy of the spheroid of NGC 3923 has a luminosity 
weighted effective radius of 3R e . 

Such a situation is very similar to that found for the 
isolated elliptical NGC 821 by Proctor et al. (2005) using 
longslit spectroscopy. In this case Proctor et al. (2005) found 
that NGC 821 displayed a strong gradient in both metallicity 
and age, in the sense that metallicity decreases with radius 
and age increases with radius. As in NGC 3923 the galaxy 
displays a constant value of [a/Fe] close to 0.3. NGC 3923 
and NGC 821 display similar metallicity gradients of ~-0.7 
dex/AlogR but the normalisation for NGC 3923 is shifted 
to higher metallicities than in NGC 821. As discussed in 
the case of NGC 821 (Proctor et al. 2005), this high value 
of the metallicity gradient tends to rule out merger models 
for the formation of NGC 3923, as merger models predict 
lower values for the metallicity gradient due to the violent 
mixing that takes place during mergers tending to wash out 
gradients (Kobayashi 2004). 

Another scenario considered by Proctor et al. (2005) to 
explain the formation of NGC 821 is where a minor merger 
event funnels gas to the inner region of the galaxy pro- 
ducing a minor star formation event and adding a layer of 
"frosting" to the majority older population already present. 
This star formation event presumably did not produce a 
significant number of star clusters stable enough or at large 
enough galactocentric radii for them to survive and be de- 
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Figure 9. Ages, metallicities and [ct/Fe] ratios measured for the NGC 3923 GCs (filled circles have S/N per A>30, unfilled have S/N 
< 30 per A), GC data binned by colour (coloured circles), NGC 3923 central regions from Thomas et al. (2005) (filled square), Denicolo 
et al. (2005) (filled triangle) and NGC 3923 spheroid (unfilled star). 



tectable at the present time. This scenario for the forma- 
tion of NGC 821, can also be applied to the formation of 
NGC 3923, where it can naturally explain the formation of 
the shell structures. However the same problem that caused 
Proctor et al. (2005) to reject this explanation in the case 
of NGC 821 still applies to our NGC 3923 data, namely 
that the gas forming the new stars would presumably be 
of different metallicity and a-element abundance than the 
galaxy at these radii- it would then have to form stars over 
the correct timescale, with the appropriate enrichment such 
that they appeared to fall exactly on the trends of metal- 
licity and a-element abundance that the rest of the galaxy 
exhibits. This appears to require an unlikely amount of fine 
tuning, however recent work by Serra & Trager (2007) shows 
that when examining a composite stellar population using 
SSP models, it is the young component which dominates 
the age determined for the composite. They further find 
that the chemical composition (metallicity and [a/Fe]) of 
the best fit SSP model closely follows that of the older pop- 
ulation which dominates the mass. This explanation would 
also not rule out the young age of around 1 Gyr determined 
by Sikkema et al. (2006) for the shell features of NGC 3923, 



as any 1 Gyr stellar population present in the inner regions 
could, with the correct choice of mass fraction, appear to be 
around 3 Gyr, when combined with a majority older stellar 
population. 

A final possible formation scenario considered by Proc- 
tor et al. (2005) to explain the formation of NGC 821 is that 
the gas for the recent starburst is fuelled by gas from the 
galaxy itself, perhaps through gas cooling smoothly onto an 
inner disc. High spatial resolution IFU observations of the 
inner regions of NGC 3923 would help to search for the sig- 
natures of such a structure. Serra & Trager (2007) note that 
if high enough S/N spectroscopy is available it is possible to 
disentangle the different stellar populations present through 
examination of the discrepancy between the SSP parameters 
determined using different Balmer-line indices. This result 
therefore provides a possible route to determining the rele- 
vance of these models for the formation of NGC 3923. 



6 CONCLUSIONS 

We have obtained low-resolution spectra of 37 GC candi- 
dates around the shell elliptical NGC 3923 using the CMOS 
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on Gemini-South. By examination of the radial velocities 
of the targets we find that 29 are GCs associated with the 
target galaxy. Making use of a new technique of extracting 
integrated spectra at extremely low surface brightness levels 
we have determined kinematic properties for the diffuse light 
of NGC 3923 out to around 200 arcsecs on both major and 
minor axes, corresponding to half light radii in the range 
2-4R e . We find evidence for low amplitude rotation along 
the major axis of around 30kms _1 . The detailed kinematic 
properties of NGC 3923 (including its observed velocity dis- 
persion in this range) and its GC system will be studied in 
more detail in a forthcoming paper. 

We have measured the ages, metallicities and [a/Fe] 
ratios of the 29 GCs and compared them to those found 
found for the diffuse light of NGC 3923 at the same pro- 
jected radii. We find all of our GCs to have old ages (>10 
Gyr), with metallicities running from [Z/H] = -1.8 to +0.35 
and a-element abundance ratios generally consistent with a 
constant value of [a/Fe] = 0.3. We see no evidence in our 
data for the existence of a younger population of GCs asso- 
ciated with the ~1 Gyr old merger event proposed to have 
produced the shell system of NGC 3923. The diffuse light of 
NGC 3923 is observed to have properties indistinguishable 
from those of the more metal rich (and redder) GCs. 

Our results provide support to the theory that the 
spheroids and GC systems of galaxies are produced dur- 
ing the same star formation events, and that the study of 
the stellar populations of GC systems can provide important 
insights into the stellar populations of galaxies at galacto- 
centric radii which at present are difficult to study using 
direct measurements of their integrated light. 
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Table 2. Measured parameters for spcctroscopically examined objects, both confirmed GCs around NGC 3923 as well as fore and 
background contaminants. In addition we provide the binned GC data and the measured indices for the spheroid of NGC 3923. For the 
GCs: ID (and number of objects per bin for the binned data in parentheses), coordinates, g magnitude, g-r and g-i colours are from our 
GMOS imaging and are instrumental magnitudes based on the standard Gcmini/GMOS zeropoints. Heliocentric velocities are from the 
spectra presented in this work, as are the measured line strength indices. For the binned GC data the mean g-i colours of the binned 
data is provided, as is the error weighted mean line strength indices for each bin, the errors for the line strengths are the 1<t scatter 
in the individual line strengths in each bin. For the spheroid data the line strengths are provided as well as errors determined from a 
combination of the statistical uncertainty on the spectrum as well as the uncertainty in the velocity determination. Data are presented 
here only for the main indices; the complete table of measured indices is available with the electronic version of this article. 
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Table 3. Derived stellar population parameters for the globular clusters and spheroid of NGC 3923. Age, [Fe/H], [o/Fe] and [Z/H] 
derived using the \ 2 minimization method described in Section 4.2. 



ID Age [Fe/H] [a/Fe] [Z/H] S/N Comments 

(Gyr) (dex) (dex) (dex) 



228 


10.0 ±1° 


-2.25±0.37 


0.48±0.40 


-1.80±0.17 


24.4 




253 


15.0 +f '! 

— o.o 


-2.28±0.34 


0.80±0.30 


-1.53±0.29 


22.5 




197 


6-0 t'i 


-0.49±0.17 


0.44±0.10 


-0.08±0.13 


29.4 




279 


15.0 +H 

— o.y 


-0.57±0.19 


0.42±0.13 


-0.18±0.16 


22.2 




112 


7.9 +\\ 


-0.88±0.13 


0.30± 0.08 


-0.60±0.15 


39.4 




394 


11.9 +H 


-1.31±0.27 


0.38±0.41 


-0.95±0.28 


15.6 


Fit unstable. 


106 


n.9 +\- 9 7 


-1.19±0.12 


0.15±0.12 


-1.05±0.08 


44.7 




64 


11.9 +_i° 7 


-1.37±0.10 


0.21±0.10 


-1.18±0.06 


56.1 




86 


15.0 +n 


-1.15±0.15 


0.34±0.11 


-0.83±0.12 


47.7 




311 


6-7 +H 


-0.72±0.37 


0.18±0.22 


-0.55±0.27 


20.0 




333 


15.0 +f 8 ° 


-0.27±0.26 


0.21±0.16 


-0.08±0.30 


17.8 


No Balmer lines in fit. 


247 


14.1 +fl 

— o.o 


-1.35±0.28 


0.32±0.24 


-1.05±0.21 


22.9 




99 


12.6 +H 


-0.57±0.11 


0.15±0.10 


-0.43±0.08 


43.4 




322 


11.2 +f 2 


-1.25±0.26 


0.24±0.30 


-1.03±0.20 


19.3 




498 


15.0 

— f .O 


0.13 ±0.25 


-0.30±0.12 


-0.15±0.29 


12.4 


Extremely poor fit 


221 


10.6 + " 


-1.45±0.25 


0.24±0.26 


-1.23±0.11 


24.9 




110 


15.0 +_H 


-0.76±0.15 


0.30±0.10 


-0.48±0.11 


40.6 




232 


H-9 ill 


-0.78±0.18 


0.24±0.15 


-0.55±0.16 


22.9 




225 


1 +4.8 
119 -3.4 


-0.85±0.24 


0.32±0.36 


-0.55±0.24 


25.7 




360 


i n /" +5.9 

10.6 +% .g 


-0.77±0.24 


0.53±0.20 


-0.28±0.14 


17.0 




332 


i k n +8.6 
15 -° -5.5 


0.15±0.16 


0.03 ±0.12 


0.18±0.18 


18.7 


Fit unstable. 


108 


14.1 +2-2 


-0.45±0.07 


0.24±0.05 


-0.23±0.05 


41.9 




450 


io.o tt 2 8 


-0.61±0.28 


0.38±0.19 


-0.25±0.22 


13.2 




513 


+ 12 

13.3 ^ 3 U 


-0.58±0.28 


0.40±0.22 


-0.20±0.23 


11.4 




492 


8.9 +H 

— 4.2 


-0.23±0.24 


0.62±0.18 


0.35±0.25 
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